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Abstract 

Using the exact Lemaitre-Bondi-Tolman solution with a non-vanishing cosmological constant A, 
we investigate how the presence of a local spherically-symmetric inhomogeneity can affect apparent 
cosmological observables, such as the deceleration parameter or the effective equation of state of 
dark energy (DE), derived from the luminosity distance under the assumption that the real space- 
time is exactly homogeneous and isotropic. The presence of a local underdensity is found to produce 
apparent phantom behavior of DE, while a locally overdense region leads to apparent quintessence 
behavior. Our study shows how observations in an inhomogeneous ACDM universe with initial 
conditions compatible with the inflationary beginning, if interpreted under the wrong assumption 
of homogeneity, can lead to the wrong conclusion about the presence of "fake" evolving dark energy 
instead of A. 
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I. INTRODUCTION 

High redshift luminosity distance measurements 1-6| and the WMAP measurement 0,0] 
of cosmic microwave background (CMB) interpreted in the context of standard FLRW cos- 
mological models have strongly disfavored a matter dominated universe, and strongly sup- 
ported a dominant dark energy component, giving rise to a positive cosmological accelera- 
tion. As an alternative to dark energy, it has been proposed jo], |lo| that we may be at the 
center of an inhomogeneous isotropic universe described by a Lemaitre-Tolman-Bondi (LTB) 
solution of Einstein's field equations, where spatial averaging over one expanding and one 
contracting region is producing a positive averaged acceleration an, but is has been shown 



how in general this procedure can lead to formal definition of unobservable quantities 
Another more general approach to niap luminosity distance as a function of redshift Dl{z) 
to LTB models has been proposed 12|, showing that an inversion method can be applied 
successfully to reproduce the observed Dl{z). Interesting analysis of observational data in 
inhomogeneous models without dark energy is given for example in [isl [l^ . 

The main point is that the luminosity distance is in general sensitive to the geometry of the 
space through which photons are propagating along null geodesies, and therefore arranging 
appropriately the geometry of a given cosmological model it is possible to reproduce a given 
Dl{z). For FLRW models this corresponds to the determination of I^a and Qm and for LTB 
models it allows to determine the functions E(r), M{r),tf,{r). 

The proposal to use galaxy number counts 15[ to distinguish between LTB models with- 
out cosmological constant and ACDM has been recently studied both analytically Ig], 



and numerically 



18| , showing how LTB models with a weak central singularity could in 



principle not be distinguished even using both the redshift spherical shell energy mn{z) and 
DLiz). 

In this paper we will take a different approach to the study of inhomogeneities, and in- 
stead of proposing them as an alternative to dark energy, we will consider they effects in 
presence of a cosmological constant, studying LTB solutions which are only locally inhomo- 
geneous whose geometry is very closed to a ACDM model, showing how even small amplitude 
inhomogeneities compatible with the amplitude of the curvature perturbation after inflation 
can lead to important effects. 

Since the amplitude of the inhomogeneities we consider is very small, corresponding 
to a few percent in terms of the density contrast, the luminosity distance Dl{z) is not 
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significantly affected, but the apparent cosmological observables derived from the Dl{z) 
under the assumption of homogeneity are significantly affected because they are sensitive to 
its derivatives. 

We consider different types of models of local inhomogeneities and obtain that local un- 
derdensity gives rise to apparent phantom behavior, while local overdensity leads to apparent 
quintessence behavior. 

Our study shows how observations of a quasi-ACDM universe with compensated large 
scale inhomogeneities compatible with inflation predictions for curvature perturbations, if 
interpreted under the '"wrong 'assumption of homogeneity, can lead to the wrong conclusion 
of the presence of "fake" evolving dark energy, while only the cosmological constant is present 
in reality. 



II. DERIVING THE EXACT LTB SOLUTION WITH A COSMOLOGICAL CON- 
STANT 



The LTB solution can be written as 19l-l2l| as 



^.^ = -rft^ + i^;^ + i?w, (1) 

where -R is a function of the time coordinate t and the radial coordinate r, E{r) is an 
arbitrary function of r, and R r = drR{t,r). 

The Einstein equations with dust and a cosmological constant give 

fR\' _ 2E{r) , 2M(r) , A 
" i?2 + i?3 +3' 

= -TrTB-^ (3) 

It ri^r 

with M{r) being an arbitrary function of r, i? = dtR{t, r) and c = SnG = 1 is assumed 
throughout the paper. Since Eq. ([2]) contains partial derivatives respect to time only, 
its general solution can be obtained from the FLRW equivalent solution by making every 
constant in the latter one an arbitrary function of r. 

The general analytical solution for a FLRW model with dust and cosmological constant 
was obtained by Edwards {22 1 in terms of elliptic functions. By an appropriate choice 
of variables and coordinates, we may extend it to the LTB case thanks to the spherical 
symmetry of both LTB and FLRW models, and to the fact that dust follows geodesies 
without being affected by adjacent regions. 



The Friedmann equation for the scale factor apit) of a pressureless FLRW universe with 
a cosmological constant has the form: 



-) =-4 + ^ + i (4) 

■.dp/ dp OCbp o 

It is convenient to introduce the conformal time rj such that drj = dt/ap, in terms of which 
the solution satisfying the initial Big-Bang condition ap{0) = can then be expressed as 

where (j){x\ Qz) is the Weierstrass elliptic function satisfying the differential equation, 
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dx 



- 92<P - 93 ■ (6) 



and we have explicitly introduced the length scale L to make the equations dimensionally 
consistent. We note that in 22[, the curvature parameter k is normalized to /c = ±1, hence 
/c^ = 1. However, for our purpose below we present the solution without normalizing k"^ to 
unity. 

We can now use this solution to construct a general solution of the partial differential 
equation ([2]). First, we introduce a new coordinate r] = ri(t,r) and a variable a by 

(^] =L^^ (7) 
[dtj^ R~ a' ^ > 

and introduce new functions by 

, , 6M(r) , , , 2Eir) 
Po(r) = ^, k{r)^—^. (8) 



1^1 =-%y + 4^a + ^a\ (9) 



Then Eq. ([2]) becomes 

2 

dr]) '^'''^ ' 3 " ' 3 
where a is now regarded as a function of r] and r; a = a{ri,r). It should be noted that 
the coordinate rj, which is a generalization of the conformal time in a homogeneous FLRW 
universe, has been only implicitly defined by Eq. ([71). The actual relation between t and rj 
can be obtained by integration t = J adrj once a(?7, r) is known. 

Inspired by the construction of the solution for the FLRW case, we can now set 

3^( jL; 93) + 
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which leads to the Weierstrass differential equation ([6]) for the choice of the parameters given 
by 

a = Po(r)L^ g2 = hirfL\ = ^ {2k{rf - Apo{rf) . (11) 
We finally get 

a{r],r] = -, r . (12) 

3<p{§;;g,ir),g,ir))+k{r)L^ 

where we have introduced the length L for dimensional consistency. In the numerical calcu- 
lations we will set L = Hq^ = 1. 



III. GEODESIC EQUATIONS 



We adopt the same method developed in 23|] to solve the null geodesic equation written 



in terms of the coordinates (?7,r). Instead of integrating differential equations numerically, 
we perform a local expansion of the solution around z = corresponding to the point 
(to,0), or equivalently (?7o,0), where tg = ^(^70)0). The change of variables from (t,r) to 
(?7,r) permits us to have r.h.s. of all equations in a fully analytical form, in contrast to 
previous considerations of this problem which require a numerical calculation of R{t, r) from 
the Einstein equation Thus, this formulation is particularly suitable for derivation of 
analytical results. 

The luminosity distance for a central observer in the LTB space-time as a function of the 
redshift z is expressed as 

Dl{z) = (1 + zfR {t{z),r{z)) = (1 + zfr{^)a {v{z),r{z)) , (13) 

where (t{z),r{z)^ or (^{ri{z),r{z)^ is the solution of the radial geodesic equation as a function 
of z. The past-directed radial null geodesies is given by 

dr . Il + 2Eir) ' 



In terms of z, Eq. ( 1141) takes the form 2J|: 



dr ^l + 2E{r{z))_ 

(15) 



dz {1 + z)RrHz),t{z)] 
dt R^r[r{z),t{z)] 



dz {l + z)R^r[r{z),t{z)] 
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The inconvenience of using the {t, r) coordinates is that there is no exact analytical solution 
for R{t, r). So the r.h.s. of Eqs. (1151) cannot be evaluated analytically, but we are required to 
find a numerical solution for R first 251], and then to integrate numerically the differential 
equations, which is quite an inconvenient and cumbersome procedure. Alternatively one 
may derive a local expansion of R{t, r) around [to, 0), corresponding to the central observer, 
and use it Eqs. (fT5l) . But one would need to expand it to a higher order in z in order to 
maintain the accuracy at high redshifts. 

For this reason, it is useful for many numerical and analytical applications to write the 
geodesic equations in terms of the coordinates (r/, r). It follows from the definition ([7]) that 



t{r],r) = tb{r) + / a{r]',r)dr] 



(16) 



hence, 



dt = a{ri, r)drj + 



^ da(r]' ,r] 
dr 



dr] + t^(r) I dr 



(17) 



Partial derivatives are transformed using the relations: 



dt 



dri^ 



Then Eqs. (fT5|) take the form: 

drj 

dz 
dr 
dz 



where 

F(?7,r) = 



R. 



2E(r) 



dr 



dr 



'df 
dr 



'd_' 
9v, 



drt{r],r) + F{ri,r) 
' (l + z)9,F(r/,r) 

1 + z)dnF{r],r) 



= p{r],r) 
q{ri,r) , 



(18) 



(19) 
(20) 



k{r)r'^ 



dr{a{ri, r)r) — a ^dj^{a{ri, r)r) drtij], r) 



(21) 



It is important that the functions p, g, F have explicit analytical forms. 



IV. APPARENT COSMOLOGICAL OBSERVABLES AND "FAKE" DARK EN- 
ERGY 

In this section we will briefly introduce the concept of apparent observables, which are 
deduced from observations assuming a fiat ACDM. 
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In a flat FLRW model for a given observed Dl{z) we have the relations, 



d (Dl{z) 
dz 



Q 



q^P^iz) 



app( 



l + Z 

rfln(g"PP(2)) 



w 



app, 
DE\ 



dln(l + z) 
(2(1 + 2) /3) d\nH^PP/dz 



q^P^{DL{z)) 
1 



(22) 
(23) 
(24) 
(25) 



i-{H^/H-ppfn^^{i + zf ■ 

These are valid under the assumption of flatness and homogeneity. Apparent observables 
are those given above with the luminosity distance Dl{z) obtained for a central observer in 
a LTB model. If the Universe is really inhomogeneous the apparent observables above will 
include the errors due to ignoring the inhomogeneity, which could for example be mistaken 
as dark energy with a redshift dependent equation of state, that is, we may be fooled by 
'fake' dark energy. 



In this paper in addition to the above, we will also consider Om{z) 26|, l27|], a diagnostic 
which can be used to distinguish ACDM from other DE models without directly involving 
the cosmic equation of state. 



Om: 



app. 



X) = 



h\x) 



X^ — 1 



X 



l + Z , h{x) = H''PP{x)/Ho . 



(26) 



V. EFFECTS OF LOCAL INHOMOGENEITIES ON APPARENT OBSERVABLES 

In order to make a connection between the LTB model and a universe with primordial 
curvature perturbations from inflation, we introduce the following metric which describes a 
spherically symmetric space-time after inflation at scales much exceeding the Hubble one: 

ds^ = -dt" + al{t)e^^^''\dr^ + r^c/fi^) . (27) 

According to the inflationary scenario, C,{r) is just a local, space-dependent number of e-folds 



constitutes the basis of the so-called 5N formalism was flrst obtained in 



in 


28| 


29 


30 



281 in case of a 



30[. Matching 



this metric to the LTB metric ([T]) with tf, = in the limit t — )■ and putting R = aF{t)e^r, 
we flnd the exact relation: 

1 + 2E(r) = [l + rC'(r)]l (28) 



In the linear approximation, this reduces to 

C'(r) 

k{r) = -2-^ . (29) 

In particular, k{0) = —2("{0). Thus, a local maximum of (, i.e. a point where inflation 
lasted longer in terms of In a, corresponds to A; > 0, i.e. to a local matter overdensity after 
inflation. Note that the LTB metric ([1]) or the relation ( !28l) is invariant under the change of 
the radial coordinate r ^ f = g{r), where g{r) is an arbitrary function of r as long as it is 
monotonic in r. We fix this degree of freedom by setting po =constant. 

Motivated by observations we consider the curvature perturbation ({r) of the amplitude 
~ 5 X 10^^. Specifically we study the four different types of inhomogeneities. 

Type r : k{r) = —[l{r] 



Type ir : C(^) = ^IKr 
Type 11+ : ({r) = A[l{r] 

where the function l{r) is defined as 

'Ho{r - ro 



l'{0)re' 






= 0.1, A = 


10-\ A = 0.02, (30) 


/'(0)re" 




HqTq 


= 0.1, A = 


-10-\ A = 0.02, (31) 


r(0)re- 


-r/r-oy 


Horo 


= 0.2, A = 


5 X 10-^ A = 0.05, (32) 


r(0)re" 


-r/roy 




= 0.2, A = 


-5 X 10-^ A = 0.05,(33) 



/(r) 



tanh 



(34) 



A ; 

In all cases, the cosmological constant is assumed to be the same as the one implied by the 
best fit ACDM model corresponding to I^a = 0.7, and Hq = H°-pp{z = 0) is adjusted to the 
observed Hubble constant. 

We have chosen these four different inhomogeneity profiles because they correspond to 
compensated inhomogeneities, in the sense that they asymptotically approach a flat, homo- 
geneous ACDM model and correspond, respectively, to 

• Type 1+ : central overdense region. 

• Type I~ : central underdense region. 

• Type 11+ : intermediate overdense region. 

• Type 11^ : intermediate underdense region. 

An important feature of type I models is that they are by construction regular at the center, 
since the linear term in the series expansion at the center is removed, avoiding the cusp 
singularity which would otherwise arise. 
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As seen from Figs. [5] and [8] for type models, the presence of a local underdensity 
gives rise to apparent phantom behavior, while that of a local overdense region to apparent 
quintessence behavior. For a shell-like underdensity or overdensity region. Figs. [H] and [IT] 
for type 11''= also indicate that an underdensity shell mimics phantom behavior, while an 
overdensity shell mimics quintessence behavior. 

Since all the models considered here describe compensated inhomogeneities, they show 
both phantom and quiescence behaviors, but we can see that the sign of the variation of 
w'de with respect to Wtme = — 1 is roughly the same as the sign of the density contrast. 

Our results give a semi-realistic example of inhomogeneities compatible with observations 
which, if interpreted in the framework of a flat and inhomogeneous spacetime, can lead to 
the wrong conclusion that there exists dark energy with an evolving equation of state, while 
in reality there is only a cosmological constant. 

VI. ANALYTICAL DERIVATION OF Qq^^ AND Wq^^ 

The calculation of the central value of the apparent cosmological observables requires to 
find an analytical expression for the right hand side of the geodesies equations. For this 
purpose we expand the relevant functions as 




(35) 
(36) 
(37) 




to get 



app 




(38) 



app 



a ri 



(39) 




The solution of the geodesies equations leads to 








(40) 



(41) 



(42) 



a' 
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We can then also expand the energy density around the center, 

n(r^ = ^° ^ rpp [Ai (r/)a,^(r/, 0) - 4a(r/, 0)a,,(r/, 0)] / 2\ /4o^ 
' a(r/,0)3^ a(r/,0)5 + j . 

In order to avoid a central singularity the term linear in r should vanish, and from this 
condition we obtain 

= 4afaO)a,faO)_ 

Using this equation we finally get 

go = ^ r ^' (^^) 



= — ^ r — / / o\ (46) 



a4 - a2 



where the right-hand sides are evaluated at (?7, r) = (?7o,0). It may be worth mentioning 
that the formulas derived so far are general in the sense they do not depend on the explicit 
form of the solution. 
We can also define 

_ a(to, 0)a(to, 0) . , 

where the derivative respect to t is denoted with a dot, and is calculated using the analytical 
solution a(?7, r) and the derivative respect to rj, 

a = d,a=^. (48) 

It is interesting to observe that because of the regularity condition we have imposed at the 
center, we have 

qr = Qo , (49) 

which can be verified using the analytical solution both for the case of vanishing and non 
vanishing cosmo logical constant. 

Using again the condition ki = we can now substitute the analytical solution to get the 
final results expressed directly in terms of observables: 

Qo - 2 M ~ ^ ^ ' W^J 

^app ^ -5Ql,CM + ^MCss{0? + Css{0f ^5^^ 

3(1 — ilM)^M 
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where we have used 



°° - C..(O)L2 + 30 ' ^^^^ 
^0 = , (53) 

^- yC..(0)^L^ -1 (16C.(0)^ + Ap^) L^) , (54) 

(55) 



( X- y Crr(0)2L^ (l6C„.(0)3 + Apt) L' 



^ 2L 



A = 3{1 -Qm)H^, (56) 

po = SalnMH^ . (58) 
As expected the above formulae reduce to the AC DM case in the central fiat limit, 

ko =Crr(0) =0, (59) 
qr = q^^'""'' = - 1 , (60) 

= wf"^^' = -1 . (61) 

As a confirmation that large scale inhomogeneities look lik;e fak;e dark; energy we can also 
verify that the relation between Qq^^ and Wq^^ is the same as in the case of an FLRW model 
with dark energy: 

go^^^^ = l^m-l + - w^il - ^m) , (62) 
qr = ^^^^ - 1 + ^(1 - wDil - ^m) . (63) 



VII. CONCLUSIONS 



We have investigated how the presence of a local inhomogeneity could affect the apparent 
equation of state of dark energy under the "wrong" assumption of a homogeneous FLRW 
background, which is commonly used in interpreting astrophysical observations in ACDM 
models. Our calculation shows how phantom and quintessence behaviors can be produced 
for compensated underdense or overdense regions. The presence of a local underdensity 
gives rise to apparent phantom behavior, while that of a local overdense region to apparent 
quintessence behavior. 
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Our results give a semi-realistic example of inhomogeneities with the amplitude com- 
patible with inflationary predictions which, if interpreted in the framework of a flat and 
inhomogeneous spacetime, can lead to the wrong conclusion of the presence of dark energy 
with an evolving equation of state. In general, a local inhomogeneity can lead to a confusion 
between local gravitational redshift and cosmological redshift due to the expansion of the 
Universe. 

Recent analysis of observational data 3l| could support the existence of a local under- 



dense region, but which many not be of compensated type as the one we considered here. 
We will investigate in a future work what could be the constraints on the size and density 
contrast of such a void based on observational data. 
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Appendix A: Calculating the density contrast 

In the text, we have carried out all our calculations in the coordinates {rj, r) since this 
allows to take full advantage of the existence of an analytical solution. But if we are 
interested in the radial profile of a quantity on a fixed time-slice t =constant, we need 
to go back to the coordinates {t,r). Below we carry this out for the density contrast, 
6 = {p{t,r) — p{t, oo))/p(t, oo), where our LTB model is assumed to approach a fiat FLRW 
universe as r — > oo. 

We need to introduce the inverse of the function defined in eq. (flGll . i.e., we need to 
express ry as a function of {t,r), r] = v(t,r), from 

t = u(ri,r) = th(r) + / a(r] ,r)dr] , (Al) 

Jo 
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such that 

u{v{t,r),r) = t. (A2) 
The value of v{t,r) can be evaluated numerically by solving for x the equation 

u{x,r) = t. (A3) 

The function 77 = v{to, r) thus obtained is plotted for the different models in Figs. l4land[T3l 
As it can be seen, rj = v{to,r) varies substantially in the region of inhomogeneity, while it 
levels off to a constant far from the inhomogeneity. 

The energy density in the coordinates {t, r) is given by 

But since the analytical solution is given in terms of 77 we have another expression, 

^(^, ,) = ^AM. . (A5) 

' a(?7, r)^r^ [9^(0.(77, r)r) — a ^9^(a(?7, r)r)9rt] 

Then the density contrast on the hypersurface t = to is given by 

c/, N p(r/(to,r),r) -p(r/(to,oo),oo) 

o{to,r) = — ^ , A6 

where 

to = t{vo, 0) . (A7) 

The density contrast is plotted in Fig. [5] for type inhomogeneities and Fig. [H] for type 
11^ inhomogeneities. 
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FIG. 1: k{r) in units of Hq is plotted for inhomogeneity of type I and I"*", r is in units of Hq ^ . 
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FIG. 2: k{z) is plotted for inhomogeneity of type I and I"*". 
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FIG. 3: C,{r) is plotted for inhomogeneity of type I and I''". 
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FIG. 4: 77 = v{to, r) in units of is plotted for inhomogeneity of types I and I"*". 




■5('or) 



0.03 
0.02 
0.01 



0.4 



0.6 



0.8 -0.01 
-0.02 
-0.03 
-0.04 




0.4 



0.6 



0.8 



FIG. 5: 5{to,r) is plotted as a function of r for inhomogeneity of types I and I"*". 
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FIG. 6: H"'PP{z) is plotted for inhomogeneity of types I and I"*". 



17 




FIG. 7: q"'PP{z) is plotted for inhomogeneity of types I and I"*". 




FIG. 8: w°jj^{z) is plotted for inhomogeneity of types I and I"*". 




FIG. 9: Om°''^^{z) is plotted for inhomogeneity of types I and I"*". 
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FIG. 10: The same as Fig. [T]but for inhomogeneity of types II and 11"*". 
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FIG. 11: The same as Fig. [2] but for inhomogeneity of types II and II"*". 




FIG. 12: The same as Fig. [3] but for inhomogeneity of types II and II"*". 
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FIG. 13: The same as Fig. |5]but for inhomogeneity of types II and 11^. 




FIG. 14: The same as Fig. [5] but for inhomogeneity of types II and 11"^. 
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FIG. 15: The same as Fig. [6] but for inhomogeneity of types II and 11"^. 
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FIG. 16: The same as Fig. [7] but for inhomogeneity of types II and 11"*". 




FIG. 17: The same as Fig. [8] but for inhomogeneity of types II and 11^. 




FIG. 18: The same as Fig. [9] but for inhomogeneity of types II and II+. 
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FIG. 19: 5Dl{z) = [D^'^^^^z)-Dl{z)]/D^'^^^^{z), the relative difference between the luminosity 
distances, is plotted for model 11^. As it can be seen the large scale inhomogeneities considered 
have a very small effect. 
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